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INTRODUCTION 


This report describes a methodology to solve the inverse design 
problem for cascades of airfoils, assuming the flow to be two- 
dimensional and the gas inviscid, compressible and ideal. The metho- 
dology is based on the procedure described in [1] for the computation 
of transonic nozzle and plume flows, in [2] for the solution of the 
inverse problem in channels, and in [ 3 ] for the inverse cascade prob- 
lem. 


The basic idea is the following. A time-dependent computation 
is performed, where a distribution of pressure is prescribed on a 
wall, the geometry of which is unknown and has to be determined. 
Such a wall is a boundary of the flow field and it is assumed as a 
flexible and impermeable surface. If the flow tends to produce a 
pressure on the wall, different from the one which was prescribed, 
the wall itself will move, adjusting its shape, until any mismatching 
is eliminated. 

The inverse problem may be formulated in different ways, both 
for airfoils and cascades. Such differences stem from a choice of 
possible prescriptions, such as the shape of the camber line, the 
thickness and/or lift distribution along the camber line, the pres- 
sure distribution on a single side of the profile, and also on the 
boundary condition assumed in front of it and behind it [ 3 ]. 

In Section 1 we report the methodology for the design of a cas- 
cade of airfoils, under the following upstream and downstream boun- 
dary conditions. Total pressure, total temperature and flow angle 
are imposed upstream. The static pressure is prescribed downstream. 
We describe the numerical technique, showing in detail how moving 
boundaries are treated for different formulations of the inverse 
problem. Numerical examples are presented for wholly subsonic flows 
and also for the design of shockless blades in the transonic regime. 
The section ends with an attempt to predict a specific shockless cas- 
cade designed by Hobson [4], by prescribing the thickness and lift 
distribution along the camber line. The method, as described in Sec- 
tion 1, fails to produce the expected solution. 

An analysis of the causes of such a failure is developed in Sec- 
tion 2. Global considerations, relating upstream and downstream 
boundary conditions with the prescribed distribution of flow proper- 
ties on the profile show that some geometries (such as the Hobson 
cascade) cannot be predicted by the methodology, unless different 
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conditions are imposed at the downstream boundary. 


In Section 3 the actual two-dimensional problem is reformulated 
in a one-dimensional approximation, introducing a suitable force 
field. Simple numerical experiments are performed to validate the 
conclusions of Section 2. A different treatment of downstream boun- 
dary conditions is proposed, as a possible way for overcoming the 
difficulties encountered with the Hobson cascade. Actual applica- 
tions of the concept to the two-dimensional problem will be discussed 
in a forthcoming Report. 

In Section 4 we present an extension of the method to the design 
of high-staggered cascades. When the stagger is high, the computa- 
tional grid described in Section 1 tends to depart too much from an 
orthogonal pattern, and accuracy can be impaired. A modified H-grid 
is then introduced, together with proper definitions of upstream and 
downstream boundary conditions. The Section ends with a presentation 
of numerical results. 


1. GENERAL FORMULATION OF THE PROBLEM 

For a cascade of airfoils, the inverse problem consists of find- 
ing the geometry of a cascade producing a flow of which some parame- 
ters are prescribed. There is a certain freedom in the formulation 
of the problem. For example, in addition to suitable conditions at 
infinity, one may: 

1) prescribe the distribution of thickness and pressure jump 
along the chord of a profile, and inquire for the geometry of the 
camber line, 

li) prescribe the distribution of thickness and pressure on one 
side of the profile, and again inquire for the geometry of the camber 
line, or 

ill) prescribe the pressure distribution around the profile 
and inquire for its geometry. 

Novak and Haymann-Haber [5] have given a solution to problem i), 
based on the Taylor expansion of the equations of the steady motion 
for a compressible inviscid flow. As they point out, such a pro- 
cedure may fail if stagnation points are in the flow field. We will 
begin by solving Novak's problem by a different technique. There- 
fore, we will consider profiles having sharp leading edges and we 
will assume that there are no stagnation points. Our technique, how- 
ever, allows more general problems to be solved, where the profiles 
have blunt leading edges, provided that the problem be split into two 
parts. On a small portion of the profile around the leading edge, 
the geometry of the profile is prescribed. The geometry of the rest 
of the profile will be determined to satisfy one of the sets of con- 
ditions mentioned above. Consequently, the stagnation point will 
find its proper location through a direct analysis of the leading 
edge region. The price to be paid is a more complicated mapping 
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needed to describe the leading edge region. 

The numerical technique described in the present paper follows 
the outlines described in [1] for the computation of transonic nozzle 
and plume flows, and in [2] for the solution of the inverse problem 
in channels. Briefly, a time-dependent computation is performed, in 
which the boundary conditions are imposed according to the formula- 
tion of the inverse problem, until a steady state is reached asymp- 
totically. The contours of the blades are considered as impermeable 
but perfectly deformable. An initial geometry is assumed. Since 
such a geometry is incompatible with a steady motion, consistent with 
the prescribed conditions, a transient is generated. During the 
transient, the walls of the blades change in shape, in order to 
satisfy the condition of impenetrability as well as the boundary con- 
ditions, compatible with problem i), ii), or iii) above. The solu- 
tion of the inverse problem is given by the geometry obtained asymp- 
totically. 

The method of solution, thus, consists of a numerical descrip- 
tion of a physical evolution in time. The partial differential equa- 
tions of motion are discretized as finite difference equations. 
Mathematically, the problem is a mixed initial-and-boundary-value hy- 
perbolic problem. Since the publication of [2], remarkable progress 
has been made in the numerical treatment of such problems, both in 
integration schemes and in the calculation of boundaries. Our 
present numerical scheme is consistent with the propagation of waves 
as implicitly described by the Euler equations of motion. It belongs 
to the A-family; details can be found in Refs. 6 through 9. in 
brief, the finite difference scheme is based on the use of one-sided 
differences, according to the direction of propagation of signals 
along the bicharacteristics in a transient, preserving the domain of 
dependence of each computed point. The main advantages of this tech- 
nique are its robustness and consistency with the physical phenomena. 
The accuracy is good even with coarse meshes and in regions where the 
gradients of the flow properties are strong or transonic transitions 
occur. Furthermore, the treatment of the boundaries, which is cru- 
cial in this type of problems, is simple and naturally consistent 
with the treatment of inner points, as well as with the overall phy- 
sical phenomenon to be described. 

For simplicity, the inverse problem is studied here under the 
assumption that the flow is irrotational and two-dimensional. 
Three-dimensional problems are only more complicated from a formal 
point of view. Straightforward extensions of the general technique 
to rotational flows have already been made [8]. 

Other approaches to the inverse problem, similar in some aspects 
to the one proposed here, have been developed independently in Refs. 
10 and 11. However, the formulations of the problem are different 
from the present one and the methods proposed seem to converge to a 
solution only if the initial configuration is very close to the solu- 
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tion itself. 


A very important point has to be discussed when dealing with in- 
verse problems: the well-posedness of the problem itself. Problem 
iii) is discussed in the case of a single airfoil in [12], where it 
is shown that the design data cannot be prescribed with complete 
freedom. In fact, they must satisfy some constraints which are dic- 
tated by the consistency of the data with the flow conditions at in- 
finity and by the requirement that the contour of the airfoil must be 
closed. The number of constraints depends on the way the inverse 
problem is formulated. This matter is analytically clear for in- 
compressible potential flow [ 13 , 111]. Unfortunately, there is not an 
exhaustive theory capable of prescribing the constraints that lead to 
a well-posed problem for compressible flow. In [12] a way is provid- 
ed to circumvent such difficulties, but the suggestion given there 
cannot be applied to the present cases, due to the peculiar formula- 
tion of the inverse problems. We address the subject with a dif- 
ferent attitude. We use a "physical" time-dependent technique. If 
the design data which we impose violate the constraints needed by the 
"steady" well-posed problem, we expect the computation never to be- 
come steady. Moreover, the well-posedness of a time-dependent formu- 
lation also depends on the upstream and downstream boundary condi- 
tions, as discussed in [15]. In fact, instabilities can be generated 
by upstream and/or downstream conditions that, together with the con- 
ditions prescribed along the blade surfaces, do not allow the flow to 
get stabilized but, on the contrary, amplify its unsteadiness until 
the computation blows up. In particular, one of the main goals of 
the present paper is to show how, having prescribed a certain down- 
stream pressure, there are two possible solutions to problem i), only 
one of which is stable from the viewpoint of a time-dependent tech- 
nique. 

Formulation ii) also presents an ambiguous feature which should 
be discussed. When the pressure distribution is prescribed on one 
side of the blade, together with its thickness, it cannot be said a 
priori whether that side has to play the role of a suction side or of 
a pressure side. The computation itself will select the role of the 
side on which the pressure distribution has been assigned. According 
to our numerical results, the computation selects that side as the 
pressure side. 
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1 . 1 The physical problem 


We proceed now to describe the physical process which we have 
chosen to generate the solution. Figs. 1.1 a) and 1.1 b) show typi- 
cal initial and final configurations. The flow is assumed to be con- 
fined between two consecutive blades, the arcs BC, and two parallel 
lines issuing from the leading edge and the trailing edge of the 
blades. The lines in front of the blades are denoted by AB. The 
lines behind the blades are denoted by CD. Such boundaries are as- 
sumed to be impermeable and perfectly deformable; therefore, we can 
think in terms of a flow within a channel, the geometry of which may 
change in time, although its width (measured parallel to the y-axis) 
is independent of time. Currently, the flow is computed only within 
a portion of the channel limited by two lines, parallel to the y-axis 
(AA and DD), which should be considered as permeable. Such artifi- 
cial boundaries are not essential. Using a suitable stretching of 
coordinates, the computational region could be extended to infinity 
in both directions. We prefer to make use of artificial boundaries 
for two reasons: 

1) the number of mesh points can be reduced, and 

2) different flow conditions at the inlet and outlet of the 
computational region can be simulated. Two boundary conditions are 
needed at the inlet surface. There, we prescribe the total tempera- 
ture 0® and the flow angle o, constant along the surface. One boun- 
dary condition is needed at the outlet surface. To begin with, we 
discuss the consequences of a simple choice of downstream boundary 
conditions, that is, prescribing the static pressure. In what fol- 
lows, it will be shown that such a choice may give rise to ambigui- 
ties in problem i). It will also be shown how uniqueness of the 
solution may be restored by a different choice of a downstream boun- 
dary condition. 

The inlet is considered fixed in time (AA does not change), 
whereas DD can slide upwards and downwards, maintaining a constant 
pitch. 

The flow repeats itself periodically in the y-direction; the 
period is the pitch of the cascade. 

Let y^ and y^ be the y-coordinates of the lower and upper wall 
of the channel, respectively. Therefore, 

yj = yi(x,t), Y2 = y2(x,t) 

The computational grid is defined by constant values of X and Y in 
the transformed set of coordinates: 

X = X 

Y = (y-yi)/(y2-yi) 

T = t 
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As Yj and Yj varY with t, the grid points slide along their 
X=constant lines. ConsequentlY, the length of the chord of the 
blades varies with time, and the soliditY of the cascade cannot be 
prescribed a priori, but it will result as a part of the solution. 

For simplicitY, the design data are prescribed as functions of 
the x-coordinate, instead of a curvilinear coordinate along the 
camber line; the restriction, however, can be easilY lifted. 


By normalizing the flow properties to reference values, the 
equations of motion for a two-dimensional, compressible, inviscid, 
isentropic flow are 


a. + ua + va +6a(u + v ) = 0 

t X y X y 

u. + uu + vu + aa /6 = 0 

t X Y X 


V. + uv + vv + aa /6 = 0 

t X Y y 

with 6 =(Y- 1 )/ 2 . In the (X,Y)-frame: 


a + ua + Aa„ + 6 a(Uy + u Y + v Y ) = 0 


+ uu^ + AUy + a(a^ + = 0 

'X ■ "’Y ■ “^yV 


V.J, + uVy + AVy + aa^Y,/6 = 0 


where 


A = Y^ + uY + vY 
t X Y 

\ - ■ V''U * 

y '^2 j-i ' 


Following [9], we rewrite (1.3) as follows: 

Uj - 


Vj - .5(f5-f,*2f3) 


where 


( 1 . 2 ) 


(1 .3) 


(l.it) 


(1.5) 


f . 
J 


A.R.„ 
J JX 


(J=1,2,3) 


- 6 - 


f j = - XjRj^ (j'lJ.S.e.T.S) 

and 

= u - a, X 2 = u + a, Xj = u, = k - aY^ 

X = A + aY , Xfi = A, X = A - aY , X« = A + aY 
5 y 6 7 X 8 ^ ( 1 . 6 ) 

= a /6 - u, R^ = a /6 + u, R^ = v 

R|j = a /6 - V, R^ = a /6 + v, ^5 “ ~ 2a/6, R^ = R^ , Rg = R^ 

The form (1.5) of the equations of motion is the result of three 
linear combinations of four compatibility equations on four charac- 
teristic surfaces. Each compatibility equation contains derivatives 
computed along two lines on a characteristic surface, viz, the bi- 
characteristic, along which signals propagate, and another line. The 
linear combinations above eliminate all terms not related to bi- 
characteristics or the particle path. Therefore, all express 

convections of a sort, and are correctly approximated using "upwind" 
differences. See in Fig. 1.2 a grid point and the projection of the 
characteristic conoid having its apex at the grid point. It is easy 
to see that the quantities Xj , defined above, are the contravariant 

components of the velocity of propagation along the four bicharac- 
teristics and the particle path, shown in Fig. 1.2. 



Fig. 1 .2 


1.2. Treatment of moving boundaries 

We consider three sets of design data, as follows: 

i) the distribution of thickness t(x) and of pressure jump 
between the two sides of the blades, Ap(x), 

ii) the distribution of thickness t(x) and of pressure along 
one side of the blades, p(x), 

ill) the distribution of pressure along both sides of the 
blades. 

We will describe the treatment of boundaries separately for each 
of the cases above. 
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i) If t(x) and Ap(x) are given 

Since the flow is periodic, the upper and lower boundaries of 
Figs. 1.1 can be reduced to a single boundary for a single blade, as 
in Fig. 1.3. Note that the upper part of the ABCD line in Figs. 1.1 
is the lower boundary in Fig. I.3, and vice versa. 

The arcs, AB and CD are deformable and impermeable interfaces, 
across which the pressure is continuous but the tangential velocity 
component may be discontinuous. In formulating the boundary condi- 
tions, the whole ABCD arc can be treated homogeneously. The inter- 
faces can be considered as surfaces of blades for which a vanishing 
thickness and a vanishing pressure jump are prescribed. With such 
convention in mind, we proceed to describe the technique for any 
blade surface. 



Fig. 1.11 


In Fig. 1.1| we show two grid points on two different sides of 
the blade, at the same abscissa. The velocity vector is decomposed 
along the tangent and the normal to the blade at each point. Since 
the blade is impermeable, the two flow velocities and the blade velo- 
city must have the same normal component. From Fig. 1.i| we see, 
thus, that 


^1t ^l^^r ^2t (1-7) 

where N^^cos <j)i , N2=cos <()2 • In addition, the thickness is constant 
in time; therefore. 


^1t ^2t 
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The pressure jump is also constant in time; consequently, for an 
isentropic flow, 

Where n=(Y-1)/('»’+1). Equations (1 .8) and (1.9) are the required 
boundary conditions. 


In the spirit of [8] and [9], the computation at the boundary is 
performed as follows. A certain number of cannot be evaluated 

at boundary points since they express signals propagating towards the 
interior of the flow field. Such terms are considered as unknowns in 
(1.5) and must be evaluated using the boundary conditions. At the 
two boundary points shown in Fig. 1.4, there are four unknowns of 
this kind: 

on the lower side, f^^ and 


on the upper side, f_.. and 

bu 

^8U yix>0' 


^8L ^2x<° 


^7U ^1x’° 


where the subscripts L and U refer to the lower side and the upper 
side of the boundary, respectively. We dispose of two boundary con- 
ditions, Eqs, (1.8) and (1.9). Therefore, we need two additional 
equations. Such a difficulty is produced by the non-orthogonality of 
the grid [16] and it is not related to the discretization process. 
Any arbitrary fixing of the difficulty based on some formal manipula- 
tion of the finite difference equations (for example, approximation 
of certain derivatives from the "wrong" side, or extrapolations of 
values to fictitious points outside the boundary) has a very high 
probability of t>eing incorrect, thus acting as a source of errors at 
every computational step. The missing equations must, indeed, 
describe some physical fact. Here, it will suffice to observe that 
the momentum equation along the boundary does not require Y- 
derivatives [9] to determine u^^ and hence 


^1x^^5U“^4U^ ^8U " ^7U " ° 

^2x^^5L”^4L^ ^8L ■ ^7L ° ° 


( 1 . 10 ) 


These equations are identically satisfied if all finite differences 
are replaced by partial derivatives. In the computational code, 
(1,10) are used to close the set of boundary conditions. Details can 
be found in C9]. Specifically, if y,jj>0 and ygx^^' boundary con- 
ditions are expressed by 


^5U ^8U^1x ^7L^2 x ^4L “ ^1 

^5U ^8U ~ ^7L " “ *^2 
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where the right-hand sides are computed using information 
interior of the flow field: 


from the 


“ ^4u“^^3U‘^^1x^^2U"^1u"^7U^'^^Vltx'^ 

3L"^2x ^ ^2L"^ 1 l"^ 8L ^ "^\^2tx 

*^2 ° " ^1U"^2U"^4u"^6U~^7U''^^1L‘'^2l''^5L‘^^6l'^^8L)^V^1 

*^3 ° ^1x^4U'^^7U 
'^4 ° ^2x^5L'^^8L 

Once the unknowns have been evaluated from (1.11), the new velocity 
of the boundary points, yi •!- “ 72^ time t+At is computed by in- 
tegrating the second derivatives, y^ = y^^^. These, in turn, are 
obtained by differentiating (1.7) and using (1.5): 


^2tt - •5'^^L-SL*2^I,-y2x»2L-f,L*‘'8L-‘'7L 


>-^V2txJ 


'"itt 

The factor, is discretized using the available values of 

the velocities, y^^. The new locations of the boundary points, y^ 

and y 2 , are obtained by integration of their velocities, v =v 

’ '^It 

Hence, the components, Vj and V 2 of the flow velocity at boundary 
points are computed as 


^1 = -2 = 

where y^^ and y^^ are evaluated using a finite difference approxima- 
tion from the new values of yj and y^. The other velocity com- 
ponents, Ui and Uj are obtained integrating their derivatives in 
time; for example. 


UlT 




Once u and v are obtained, the new values of u and v are easily 

puted from u and v and the new values of y, and . 

''lx '' 2x 


com- 


On either side of the blade, the new values of the speed of 
sound, a, are evaluated using the first of (1.5). 

A less cumbersome method consists of solving (1.11) and then 
computing u, v and a from (1.5). The blade velocities can be 
evaluated from (1.7). Numerical experiments, however, have shown 
that the more complicated method, described above, yields more accu- 
rate results. 


ii) If t(x) and p(x) on one side are given 


The problem is solved as described above, with the following 
change: The boundary condition (1.9) is replaced by the condition, 

p =0, that is, for isentropic flows, 

I U 

= 0 ( 1 . 12 ) 

if, for instance, the pressure is prescribed on the upper side of the 
blade. Consequently, the second of (1.11) must be replaced by 


5U 




(1.13) 


where 


0 = - f -f *f ~f -f 

^5 1U 2U m 6U 7U 

For the remaining portions of moving boundaries (the arcs AB and CD 
in Fig. 1.3) the calculation proceeds as in case i). 


iii) If p(x) is prescribed along both sides of the blade 


This problem is the same as the one described in [2] for design- 
ing the walls of a duct, with some modifications required by the dif- 
ferent numerical scheme used here. The boundary conditions have to 
provide that the pressure does not change in time on both sides of 
the blade. Therefore, the second of (1.11) is replaced by (1.13)#as 
in case ii) above and, similarly, the first of (1.11) is replaced by 




(I.IJI) 


Where 


^6 = 


f -f -f „f «f 

1L 2L 5L 6L 8L 


Again, the interfaces AB and CD (Fig. 1.3) are treated as in case i). 
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Fig. 1.6 


1.3. Numerical results 

To illustrate our method and to discuss the well-posedness of 
the formulation, we present three sets of results: Two stable calcu- 
lations for which t(x) and Ap(x) are prescribed, two calculations for 
which t(x) and Pi(x) are prescribed (one being stable and the other 
unstable), and an unsuccessful attempt to recompute the Hobson cas- 
cade No. 2 [4] according to formulation i). 

Fig. 1.5 shows the initial configuration and Fig. 1.6 the steady 
solution to the inverse problem for the case corresponding to 

T = .025 [1 - cos(2nx)] (0<x<1) 

and 

Ap = .1 [1 - cos(2irx)] (0<x<1 ) 

The ratio P^/p® between downstream pressure and total pressure is 

0.8, the upstream slope of the velocity vector, o is 20°, and the 
upstream total temperature 0“ is 1 . Both this case and the following 
one have been computed using ^10 intervals in X and 10 intervals in Y. 



Fig. 1 .7 
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A check on the accuracy of the computation is shown in Fig. 1.7, 
where the theoretical behavior of the integral of the y-momentum is 
compared with the numerical result. The maximum error is less than 
0 . 1 ?. 


The case of Fig. 1.8 has the same t, o, and 0° as in the preced*- 
ing case, but 

Ap = .15 [1 - cos(2ttx)] (0<x<1) 

and p^/p®=0.71 . 



The resulting cascade is supercritical but unchoked and shock- 
less. It can be seen from the isomachs of Fig. 1.8 that a supersonic 
bubble appears on the upper side of the blade, but the lower side is 
entirely subsonic. The pressure cannot be discontinuous on the sub- 
sonic side; therefore, it must be continuous on the supersonic side 
as well, since Ap is prescribed as a continuous function of x. The 
method does not imply that the whole flow field is shockless, but 
just that the blade surface is shockless. The behavior of the 
isomachs of Fig. 1.8 may suggest that an imbedded shock is present 
inside the flow field and it vanishes at the blade wall. The method, 
indeed, is capable of capturing shock waves, in the sense described 
in [7]. The mesh resolution of Fig. 1.8 is not sufficient to provide 
a conclusive interpretation of the result; anyway, the possibility of 
designing blades that are shockless just at the wall is of great 
relevance from a practical viewpoint, especially in the case of su- 
percritical compressor cascades. 

In this case, as well as in the preceding one, the theoretical 
behavior of the y-momentum has been compared with the numerical 
results and the maximum error is less than 0.3?. 

In the next two examples, p(x) is prescribed on one side of the 
blades. As mentioned above, this problem presents an ambiguity. It 
is not specified, indeed, whether the side along which the pressure 
is prescribed is the suction side or the pressure side. Our numeri- 
cal experiments show that, of the two possibilities, only the one for 
which the pressure is prescribed on the pressure side is stable. We 
do not intend to say that a steady, unstable flow may exist about a 
cascade with rigid blades. The instability detected in the numerical 
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solution is related to the problem as posed here, where the blades 
are perfectly deformable and the data are prescribed as said above. 

In the numerical experiments, the steady solution to the inverse 
problem i) has been used as the initial condition for the inverse 
problem ii). The pressure distribution on one side of the blades has 
been taken equal to the previously computed pressure distribution and 
forced to remain unchanged. 



Fig. 1.9 


The initial configuration and related Isobars are shown in Fig. 
1.9. They solve problem i) for 

T = .025 [1 ” cos(2irx)] (0<x<1 ) 

and 

Ap = - .075 [1 - cos(2t7x)] (0<x<1 ) 

with p^/p°=0.8, 0 = 0 , 0°=1 . 

Having forced the pressure on the upper side (pressure side) of 
the blades to remain unchanged, 1000 more integration steps were tak- 
en. The pattern of Fig. 1.9 remains unchanged. 

Then, a second computation is performed for which the initial 
condition, shown in Fig. 1.10, is the mirror image of Fig, 1,9. The 
pressure is forced again to remain unchanged on the upper side of the 
blade, which is now the suction side. 



Fig. 1.10 
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After a few computational steps, minor numerical perturbations 
trigger a violent transient, which does not converge to a steady 
solution. Fig. 1.11 shows the flow field after 800 integration steps 
and Fig. 1.12 after 850 steps. In Fig. 1.13 the interface-and-blade 
velocity distributions are plotted at steps 1000 and 2000. Obvious- 
ly, the flow, as computed, does not converge to a steady state. 

We interpret the last two examples as follows: In the case of 
Fig. 1.9 the results of problem ii) perfectly agree with the results 
of problem i). The design pressure distribution, as assumed from a 
previous computation, is consistent with a well-posed problem for the 
pressure side. In the last example, instead, we observe two facts. 
First, the initial configuration is unstable for problem ii) and a 
transient is generated in an attempt to find a different steady solu- 
tion with the upper side of the blades playing the role of a pressure 
side. Second, the assigned pressure distribution is inconsistent 
with a well-posed problem for the pressure side and the computation 
is unable to converge to a steady state. 

Finally, Hobson’s cascade No. 2 [4] and its steady flow field 
were adopted as initial conditions for a calculation of the first 
type, and we expected our computed results not to depart substantial- 
ly from Hobson’s. The evolution of the interface and blade lower 
surface, however, turns out to be very different from what we expect- 
ed, as shown in Fig. 1.14. Clearly, the blade geometry tends to 
depart more and more from the original one. The question then ar- 
ises, Why is the third calculation stable and the latter unstable, 
despite the apparent similarity of the procedure and of all the ini- 
tial and boundary conditions? We will attempt to find an answer in 
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2. FLOW DEFLECTION AND FORCE ACTING ON A BLADE 

With reference to Fig. 2.1, let us consider a cascade with inlet 
and outlet boundaries located sufficiently far upstream and down- 
stream, so that the flow at such boundaries does not depend on y. 
Let us prescribe total temperature (0^ and flow angle (o) upstream, 

and static pressure (p^) downstream. The flow is considered isentro- 

pic. In a steady state configuration, the tangential force (F) act- 
ing on one blade is related to the upstream and downstream y~velocity 
components: 
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( 2 . 1 ) 


F = m (v^ - v^) 

where m is the mass flow through the stream channel corresponding to 
the pitch (s) of the cascade. 



On the other hand, F can be determined as a function of the exit 
angle, o^, as follows: 

a) The exit velocity, q is a known function of 0.® and p : 

e 1 e 


\1 6 


26/Y,, 

Pe ) 


b) Assume a value of o , and let E = tan o , 

e* e e 

c) The velocity components, u^ and v^ are: 


u = 
e 


e 


V = 

e 


q E 
^e e 


(Ue2)1/2 

e 


d) The mass flow can thus be evaluated: 


m = s u p , 
e ^e' 


with p = p 


1/Y 


e) p^ and the inlet velocity components are obtained by solving 


m = s 




( 1 +tan^o^) 


1/2 ^i 


where 


Pi = Pi 


1/Y 


" \| 


i (Qi” - pf'") 


and m is the same as obtained in d) . 
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f) Now the force F is computed, according to (2.1). 

From a prescribed set of values for 0^®, and p^, one can 

compute f’(Og)* This relationship is plotted in Fig. 2.2. We see 

that F vanishes for three values of o . At o = ±90° , because u 
. e e e 

and, consequently, m, vanish. In these two points, the blades are so 

deflected at the trailing edge, that the flow has no axial velocity 

component. In addition, F vanishes for that is, when the 

blade does not deflect the flow. The force is positive when 

and vice versa. 



Assume now that the inverse problem has to be solved with a set 
of boundary conditions (0^°, and p^) and a prescribed distribution 

of lift over one blade (Ap(x)). On the basis of the previous con- 
siderations, two different geometries of the blade may satisfy the 
problem, or none. The force F results as 

F = /Ap(x)dx 

with the integral carried over the x-axis between leading and trail- 
ing edge. If F>F or F<F . , no solution exists, but if 
^ ^ max min 

F . <F<F , two different cascades may be obtained, providing the 
min max & 

same force (such as A and B, for example). The same force, thus, may 
be balanced, according to (2.1), by a lower mass flow and a larger 
deflection at point A, or by a higher mass flow and a smaller deflec- 
tion at B. However, if we consider the exit boundary condition as in 
the previous Section, only the configuration described by point B can 
be reached. Let us, indeed, consider a rigid cascade providing an 
exit angle, o^, slightly smaller than (o^)^, and a steady flow 

through it, which acts on the blades with a force slightly smaller 
than F^. Let us now increase the pressure jump of this initial con- 
figuration so that the force reaches the value F^ and, at the same 

time, let the blade move to adjust itself to the new condition. 
Since F^ is larger that the initial force, the curvature of the blade 
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must increase and o decreases, instead of increasing towards (o ) . 

® e A 

Thus, the blade geometry tends to move farther and farther away from 
A. The opposite motion of the blade occurs if the initial value of 
Og is slightly larger than (o^)^, and the force, originally greater 

than F^, is decreased. In this case, however, the geometry of the 

blade will eventually reach point B. In conclusion, B represents a 
stable configuration and it is the only one which can be reached us~ 
ing the numerical procedure and the exit boundary condition, as for- 
mulated above. Do the cascade shown in Fig. 1.8 and the Hobson cas- 
cade belong to geometries of type A or B? Their curves, F( 0 g), ob- 
tained for the values of 9^®, and p^ prescribed in Section 1, are 

plotted in Fig. 2.3, where the curve labeled (1) refers to Fig. 1.8 
and the one labeled (2) refers to the Hobson cascade. It Is in- 
teresting to see that the case of Fig. 1.8 is of the B-type, whereas 
the Hobson cascade is of the A-type. It is now clear why the 
geometry of Fig. 1.8 is very stable, and why we have not been able to 
get a stable configuration for the Hobson cascade, even if the compu- 
tation was started near the desired solution. In conclusion, the 
Hobson cascade cannot be obtained, unless we modify the formulation 
of the inverse problem. 



Fig. 2.3 


3. ONE-DIMENSIONAL APPROACH 

For an easier understanding of the conclusions of Section 2, we 
will resort to a one-dimensional analysis. 

The actual cascade is composed of blades of finite thickness, spaced 
apart by a finite pitch. Let us replace it by a cascade of blades of 
vanishing thickness. Infinitely close to one another. Consequently, 
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all y-gradients vanish. Instead of prescribing the pressure jump 
distribution (Ap(x)) over a single blade, we must prescribe a force 

field per unit volume (f^) along x, acting along the y-axis. 

The equations of motion, formerly written as in (1.2), become 


a^ + ua + 6au = 0 

t X X 

a n 


(3.1) 


V, + uv + — =0 

t X p 

Here f^ represents the x-component of the force per unit volume. Be- 
cause of the assumption of inviscid flow, the total force, 

f = i + J 

is, at each point and at any time, perpendicular to the instantaneous 

y 

shape of the blade, as shown in Fig. 3*1* So, while f is prescribed 

as a constant in time, f will change, during the transient, as the 
slope (b^) of the blades varies; 


Eqs (3.1) then become: 


= - b 

X 


a^ + ua + 6au = 0 

t X X 


u^ + uu + — a 

t X 6 X 


b f-" 

X 


= 0 


(3.2) 


v^ + uv + — = 0 
t X p 


Finally, the blade moves locally with the velocity, b^ defined by: 


= - ub^ + V (3.3) 

Eqs (3.2) can be integrated in time, after having been rearranged in 
the proper form, as done in Section 1 (see Eqs 1.5). Once u and v 
are known, b^ may be evaluated and a new geometry of the blade may be 

predicted. 

We may now proceed to some numerical experiments. Let us con- 
sider a cascade with the following set of boundary conditions: 
0i®=1; 0^=20° ;p^=0,9. The curve F(o^) for this case is plotted in 

Fig. 3.2. We assume the following distribution of f^ over the blade: 
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Fig. 3.1 


= 0 


( 0 <x<x 


r = 0.1238 [1 - cos 2 t7 


x-x 


l.e. ’ 
l.e. 




(3.11) 


X. - 

t.e. l.e, 


] 




The force (per unit length along y) is: 

F = /f^dx = 0.1238 


The two possible configurations for the blade, corresponding to this 
force, are denoted by points A and B (o .= -49.28° and o „= -20°). 

SH 0D 

Their geometries can be easily found. At any x, the relation between 

y 

V and f , the conservation of mass and total enthalpy yield the fol- 
lowing three equations for the three unknowns, v, u and p: 


V 


u p = m 


V. 

1 


. 1 ^f "dx 


m o 


, 1/Y. 

(p = P ) 


(3.5) 


26/Y 

P 



+ v^) 


1 


where ra is known. From u and v, the slope of the blade is obtained 
at each abscissa, x: b^= v/u, and the blade geometry, b(x) follows, 

by numerical integration. 


Let us consider first the cascade related to point B on Fig. 
3.2. Its shape is reported in Fig. 3,3 as the curve labeled "K=0". 
We assume this geometry as the initial condition of a time-dependent 
computation, by letting the blade move and not changing the force. 
In the first integration steps the blade moves slightly and then it 
gets stabilized. The geometry at K=1000 is shown in Fig. 3 . 3 . The 
slight differences between the geometries labeled K =0 and K=1000 are 
due to the different b(x) obtained through the integration following 
( 3 . 5 ) and by the numerical time-dependent solution. The same stable 
result is reached if one starts from a geometry slightly different 
from the correct one. 
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Fig. 3.2 


We consider now the cascade related to point A. Its shape is 
given by the middle line in the set shown in Fig. 3.^. Once more, 
starting from this initial condition, we perform a time-dependent 
computation letting the blade move and keeping the force unchanged. 
As shown in Fig. 3.5, the blade tends first to move slightly down-^ 
wards, seems to get stabilized there, but at K=600 starts to move up- 
wards until it reaches another stable configuration, which is the one 
corresponding to point B. If the computation starts from the 
geometry shown by the upper line in Fig. 3-^, the behavior is similar 
to the one shown in Fig. 3.5, but the convergence to point B is even 
faster, as shown in Fig. 3*6. 

Finally, assuming the lower line of Fig. 3-^ ss the initial 
shape, we find that the blade moves as shown in Fig. 3.7 and the com- 
putation blows up. Let us point out that the behavior of the latter 
numerical experiment (Fig. 3*7) is very close to the degenerative 
process shown in Fig. 1.14 for the Hobson cascade. 

In order to succeed in obtaining a geometry of the A-type, one 
should try to make the function o^(F) single valued. All difficul- 
ties, indeed, seem to arise from the fact that different mass flows 
can provide the same force with different deflections. To achieve 
our goal, we replace the downstream boundary condition (constant 
static pressure p^) by an exit surface modeled in the spirit of Ref. 

17. As shown in Fig. 3*8, a discontinuity is located at the exit 
boundary, which simulates a guide forcing the flow to be discharged 
to the right at the prescribed static pressure, p^^, and with the 

given angle, o . The flow coming onto the discontinuity 

(u^, v^, p^) will be determined by assuming that mass flow and total 
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Fig. 3.6 
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Fig. 3.7 


pressure are the same on both sides of the discontinuity: 

p u = p u 
e e ex ex 

^1/Y _ 1/Y 

^e ^e ’ *^ex ^ex 

26/Y 6,2 2, 26/Y 6 ^ 2 , 

P +tt(u+v)=p + — (1+tan 0 ) u 

e Y e e *^ex Y ex ex 

Now we can prescribe a set of boundary conditions, 0.®, o., p , o , 

1 X 0 X 0X 

and a new relationship, ^(o^) can be found: 

a) From p and o , the mass flow is computed: 
ex ex 


m = p 


ex (1+tan^o ) 
ex 


I (Pe’"*''' 


^ex ' 


b) Using (3-5), u,v and p are computed at all x-locations; in 

particular, we obtain u , v , p (and o = arctan (v /u )) at the 

© © © © © © 

exit. 


The function, F(Og) is now monotonic, as shown in Fig. 3.9. The 
curve is limited by two points, M and N. For and it is 

no longer possible to maintain the same mass flow required by p^^ and 
If a force is prescribed, such that Fj^<F<Fj^, there is only one 
acceptable geometry which can be shown to be stable, using an argu- 
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Fig. 3.8 


ment similar to the one employed above. Our conclusions can be test- 
ed numerically. In particular, if we perform the computations with 

= o. = 20°, F=0.1238 and o = -20“ or o = -49.28“, which 

11 ex ex 

correspond to the two points A and B of Fig. 3*2, we are able to 

predict both shapes very well and both results are perfectly stable. 
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i|. DESIGN OF CASCADES WITH HIGH STAGGER 


If the projected cascade turns out to have a high stagger, the 
H-grid which we have been using so far will become very skewed, as 
shown in Fig. ^.1. Orthogonal grids are known to provide the highest 
degree of accuracy, and it is generally advisable to keep the compu- 
tational grid as close to orthogonality as possible. A grid such as 
the one shown in Fig. 4.2 is better suited for numerical work than 
the grid shown in Fig. 4.1. The upper blade is shifted with respect 
to the lower one by a number (NST) of intervals, each of a length Ax, 
in the positive x-direction. Periodicity is provided by making 
points A,B,C and D correspond to points G,I,J and K, respectively. 
Periodicity conditions are imposed on AB and GI and on CD and JK. 
The prescribed thickness and the jumps in pressure are applied to 
correspondingly shifted points (for example, M and N). The inlet 
(AF) and outlet (EK) are treated as explained in the previous Sec- 
tions. In particular, the total temperature (0^®) and the flow angle 

(o) are prescribed at the inlet and the static pressure (p^) at the 

outlet. Points on the segment FG (except G) near the inlet and 
points on DE (except D) at the outlet have no correspondence with 
other points. The segment, FG is considered as a solid wall, having 
the prescribed flow direction (o) at the inlet. The segment, DE is 
assumed as a flexible wall, which will adapt its shape in order to 
maintain a static pressure equal to p^. Point G corresponds to A, 

and point D to K. The inlet blade angle, aj’, as shown in Fig. 4.3, 
is obtained as 




t 



arccos 


(NST) Ax 
P 


once the pitch, p is prescribed. Different inlet blade angles may be 
simulated by properly adjusting three parameters; 


1) the value of (NST), which defines a discrete shifting of the 
grid points of the upper blade with respect to the corresponding 
points of the lower blade, 

2) the value of Ax, which is determined by the number of points 
over the blade of a prescribed length, and 

3) the inlet flow angle. 


The freedom in choosing such parameters allows different final 
stagger angles of the cascade to be simulated. 

A typical result of the design procedure with the grid described 
above is shown in Fig. 4.4. 
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